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HD 60532, a planetary system in a 3:1 mean motion resonance. 
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ABSTRACT 

In a recent paper it was reported a planetary system around the star HD 60532, composed by two giant planets in a possible 3:1 mean 
motion resonance, that should be confirmed within the next decade. Here we show that the analysis of the global dynamics of the 
system allows to confirm this resonance. The present best fit to data already corresponds to this resonant configuration and the system 
is stable for at least 5 Gry. The 3:1 resonance is so robust that stability is still possible for a wide variety of orbital parameters around 
the best fit solution and also if the inclination of the system orbital plane with respect to the plane of the sky is as small as 15°. 
Moreover, if the inclination is taken as a free parameter in the adjustment to the observations, we find an inclination ~ 20°, which 
corresponds to = 3.1 Mjup and = 7.4 Mjup for the planetary companions. 
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1. Introduction 



In a recent paper (I Pesort et~an l2008l hereafter DES08), two 
planetary mass companions were detected around the F type star 
HD 60532. Using the data acquired with the HARPS spectro- 
graph based on the 3.6-m ESO telescope at La Silla Observatory, 
and a two-keplerian fit to the data, the orbits of the two bod- 
ies were determined, corresponding to minimum masses of 
1.03 Mjup and 2.46 Mjup at 201-day and 604-day period, respec- 
tively. The dynamical study done in DES08 suggested a possi- 
ble 3:1 resonance, but the stability of the results seemed to be 
questionable, as small variations of the semi major axis led to 
a non resonant solution (DES08). Therefore, the authors con- 
cluded that the existence of a mean motion resonance would re- 
quire an additional decade of observations to be confirmed. 

The presence of two or more interacting planets in a system 
dramatically increases our potential ability to constrain and un- 
derstand the processes of planetary formation and evolution. The 
dynamical analysis of such systems is then very useful, first for 
constraining the system evolution history and second to deter- 
mine the system "structure" in terms of orbital content. Among 
the known multi-planet systems, a significant fraction present 
strong interactions and are trapped in mean motion resonances. 
The resonances between planets are believed to be formed after 
their inward or outward migration of the planets during the ear ly 
stages of the formation of the system (eg.^T siganis et al.ir2005l) . 

By performing a detailed dynamical analysis of the 
HD 60532 system in conformity with the radial-velocity obser- 
vations, we can confirm this system in a 3:1 mean-motion reso- 
nance and provide some constraints on the inclination of the sys- 
tem orbital plane with respect to the plane of the sky. In Section 2 
we re-analyze the observational data obtained with the HARPS 
spectrograph and the dynamical analysis of the system discussed 
in Sect. 3. Finally, conclusions are drawn in Sect. 4. 
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2. Orbital solution for the HD 60532 system 

The 147 published radial- velocity data points of HD 60532 
(DES08) were taken with the HARPS spectrograph from 
February 2006 until June 2008. A two-companion Keplerian so- 
lution is provided in DES08, but because of the strong mutual 
interactions between the two planets, these parameters present 
important variations and cannot be used in dynamical studies 
unless the initial longitudes of both planets are given. Usi ng the 
iterative Levenberg-Marquardt method ( Press et al.l Il992h . we 
thus re-fitted the complete set of radial velocities with a 3 -body 
Newtonian model, assuming co-planar motion perpendicular to 
the plane of the sky, similarly to what has been done for the 
system HD 45364 ( Correia et al. 2009). This fit yields an adjust- 
ment of = 4.41 and rms = 4.38 ms"^ slightly better than 
the two-Keplerian model from DES08. The set of orbital param- 
eters for this system (Si) is given in Table[T] 

Still assuming co-planar motion, we then release the perpen- 
dicularity constraint by including the inclination as a free pa- 
rameter. The system readjusts slightly the orbital parameters of 
the two planets and provides for the best stable fit an inclina- 
tion around 20°, although subjected to a large uncertitude. Fixing 
the inclination at this value the new fit yields an adjustment of 
^|x^ = 4.37 and rms = 4.34 ms"^. In this case the masses are 
increased by a factor 1/sin/ = 2.92, i.e., = 3.1Mjup and 
Mc = 7.4Mjup (S2, Table[I]). We also attempted to fit the data 
with a 3 -body Newtonian model for which the inclination of the 
orbital planes, as well as the node of the outer planet orbit, was 
free to vary. We were able to find a wide variety of configu- 
rations, some with very low inclination values for one or both 
planets, that slightly improved the fit. However, all of these de- 
terminations are uncertain, and since we also increased the num- 
ber of free parameters by three, we cannot say that these solu- 
tions present a real improvement with respect to the sets <Si, <S2, 
presented in Table[T] 
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Table 1. Orbital parameters for the HD 60532 system Si, ob- 
tained with a 3 -body Newtonian fit with / = 90° (top) and for the 
system system S2, with / = 20° (bottom) (M^ = 1.44 M©). 



Param. 


[unit] 


HD 60532 b 


HD 60532 c 






4.407 




[111/ j> J 


4.380 


Date 


[JD-2400000] 


54000.00 (fixed) 


V 


[km/s] 


-0.0049 + 0.0004 


P 


[day] 


201.46 + 0.13 


605.28 + 2.12 


/L 




14.48 + 0.69 


316.23 + 0.95 






0.279 + 0.006 


0.027 + 0.007 


CJ 


[deg] 


352.15 + 1.08 


136.81 + 16.34 


K 


[m/s] 


29.63 + 0.32 


46.80 + 0.40 


i 


[deg] 


90 (fixed) 


90 (fixed) 


M 


UVljupJ 


1.0484 


2.4866 


d 




0.7597 


1.5822 


Param. 


[unit] 


HD 60532 b 


HD 60532 c 






4.369 




Liii/j>j 


4.342 


Date 


[JD-2400000] 


54000.00 (fixed) 


V 


[km/s] 


-0.0055 + 0.0003 


P 


[day] 


201.83 + 0.14 


607.06 + 2.07 


A 


[deg] 


14.78 + 0.66 


317.02 + 0.93 


e 




0.278 + 0.006 


0.038 + 0.008 


CJ 


[deg] 


352.83 + 1.05 


119.49 + 9.14 


K 


[m/s] 


30.34 + 0.32 


47.84 + 0.44 


i 


[deg] 


20 (fixed) 


20(fixed) 


M 


[Mjup] 


3.1548 


7.4634 


a 


[AU] 


0.7606 


1.5854 
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Fig. 1. Variation in the resonant argument, = At - 3Ac -\- 2ojb 
(a) and in the secular argument, Aoj = ojh - ojc (b), with time. 

3. The 3:1 mean motion resonance 

As in DES08, our nominal solution Si (Tab[T]) seems to be in 
a 3:1 resonance. We integrated numerically the orbits of the 
planets over 100 kyr with the symplectic integrator SABAC4 of 
iLaskar & Robutell (l200ll) , using a step size of 0.02 years. The 



Table 2. Fundamental frequencies for the orbital solution Si 
(Tabd]). Hh and ric are the mean motions, gi and g2 are the secular 
frequencies of the periastrons, and le is the libration frequency 
of the resonant angle = Ah - 3 Ac + 2ajt,. 





Frequency 


Period 




Vyr 


yr 


rib 


653.645983 


0.550757 


He 


217.563988 


1.654686 


gi 


-0.477010 


-754.701 


82 


-0.009651 


-37303.55 


le 


9.015796 


39.9299 



Table 3. Quasi-periodic decomposition of the resonant angle 
6 = Ab - 3 Ac + 2ojh for the orbital solution Si (TabU]). We have 
6 = 180 + cos(vi t + (pi), where the amplitude and phases 

Ai, (pi, are given in degree, and the frequencies v/ in degree/year. 
All terms are identified as integer combinations of the funda- 
mental frequencies given in Table[2l which is a signature of a 
very regular motion. 





Combination 






Ai 


0/ 




He 


^1 


g2 


k 


(deg/yr) 


(deg) 


(deg) 








-1 


1 





0.4674 


63.652 


-81.114 














1 


9.0158 


38.138 


168.018 








1 


-1 


1 


8.5484 


35.869 


159.131 








2 


-2 


1 


8.0811 


25.845 


-29.755 








3 


-3 


1 


7.6137 


16.868 


141.358 








4 


-4 


1 


7.1464 


10.622 


-47.528 








5 


-5 


1 


6.6790 


6.325 


123.585 








6 


-6 


1 


6.2116 


3.509 


-65.301 








-2 


2 





0.9347 


2.697 


107.773 








-1 


1 


1 


9.4832 


2.366 


176.904 








7 


-7 


1 


5.7443 


1.872 


105.812 








3 


-3 


2 


16.6295 


1.156 


39.376 








2 


-2 


2 


17.0969 


0.874 


-131.738 








-2 


2 


1 


9.9505 


0.874 


5.791 














2 


18.0316 


0.828 


66.035 








8 


-8 


1 


5.2769 


0.964 


-83.074 








4 


-4 


2 


16.1622 


0.805 


-149.510 








1 


-1 


2 


17.5642 


0.499 


-122.851 








-3 


3 


1 


10.4179 


0.517 


-165.323 








-3 


3 





1.4021 


0.433 


-63.341 








9 


-9 


1 


4.8096 


0.482 


88.040 


1 


-1 











436.0820 


0.395 


-173.181 



frequency analysis of this orbital solution allows then to con- 
clude that it is indeed in a 3: 1 mean motion resonance, with res- 
onant argument: 6 = Ah - 3 Ac + 2ojiy. 

The fundamental frequencies of the systems are the two 
mean motions nt and ric, the two secular frequencies of the peri- 
astrons gi and g2, and the libration frequency of the resonant 
argument le (Tab©. We have, up to the precision of the de- 
termination of the frequencies 10"^), the resonant relation 
rib - 3nc + 2gi = 0. 

For the nominal solution Si , the resonant argument is in 
libration around 180°, with a libration period 27i/le = 39.93 yr, 
and an associated amplitude of about 38.138 degrees, with mul- 
tiples additional terms of significant amplitude of the form 
k(gl - g2) -h le, with k integer (Fig. [3^, Table[3]). For the com- 
plete solution, the libration amplitude can reach 123.16 degrees 
because additional periodic terms are present. In Table[3j we pro- 
vide a quasi-periodic decomposition of the resonant angle 6 in 
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Fig. 2. Stability analysis of the HD 60532 planetary system for different values of the inclination (90°, 20°, 10°). For a fixed orbital 
inclination, and initial condition of the outer and inner planet, the phase space of the system is explored by varying the semi-major 
axis Gk and eccentricity ek of the other planet, (b) and (c) respectively. The step size is 0.002 AU in semi-major axis and 0.005 in 
eccentricity. For each initial condition, the full system is integrated numerically over 10 kyr a nd, as in (Correia et al. 2005, 20Q1), 
a stability criterion is derived with the frequency analysis of the mean longitude (ILaskarll 1 99QI 1 1 993]) . The "red" zone corresponds 
to highly unstable orbits, while the "dark blue" region can be assumed to be stable on a billion-years timescale. The contour 
curves indicate the value of rms obtained for each choice of parameters. It is remarkable that in the present fit, there is perfect 
correspondence between the zone of minimal rms and the 3:1 stable resonant zone, in "dark blue". 



terms of decreasing amplitude. All the quasi-periodic terms are 
easily identified as integer combinations of the fundamental fre- 
quencies (TableO. Since the resonant angle is modulated by a 
relatively short period of about 38 years, and high harmonics of 
shorter periods, the observation of the system over a few addi- 
tional years may provide an estimate of the libration amplitude 
and thus a strong constraint on the parameters of the system. 

As in (Correia et al 2009), both periastrons precess with 
mean frequency gi that is retrograde, with a period of 754.701 
years. The two periastrons are thus locked in an antipodal state, 
and the diff'erence Aoj = ojh - oJc is in libration around 180° 
with an amplitude of 90? 7 (Fig.[3]3). As a result, the arguments 
6' = Ah -3Ac -\- 2ajc, and 0'' = At - 3Ac -\- oJb ojc librate around 
180° and 0° respectively, with the same libration frequency Iq. 



4. Stability analysis 

To analyze the stability of the nominal solution and confirm the 
presence of the 3 :1 resonance, we performed a global frequency 
analysis (lLaskarl»1993.) in the vicinity of the nominal solution 

i Fig. in, in the same way as it was achieved for the HD 202206 
Correia et al.1 l2005h and HD 45364 (Correia et al." "2009') sys- 
tems. For each planet, the system is integrated on a regular 2D 
mesh of initial conditions, with varying semi-major axis and ec- 
centricity, while the other parameters are retained at their nomi- 



nal values. The solution is integrated over 10 kyr for each initial 
condition and a stability indicator is computed to be the varia- 
tion in the measured mean motion over two consecutive 5 kyr 
intervals of time. The results, for diff'erent values of the inclina- 
tion of the system (/ = 90°, 20°, 10°), are reported in color in 
Fig.m where "red" represents the strongly chaotic trajectories, 
and "dark blue" the extremely stable ones. In all these plots, the 
3:1 resonant island is clearly highlighted, and appears to be very 
stable until the inclination reaches / = 10°, for which most of the 
resonant island is destabilized. Indeed, a numerical integration of 
the fitted solution for 10 degrees inclination led to a disruption 
of the system in 584 Myr. 

It is remarkable that, as for HD 45364 (ICorreia et al.ll2009l) . 
there is perfect coincidence between the stable 3:1 resonant is- 
lands, and curves of minimal rms obtained by comparison with 
the observations. Moreover, for our preferred solutions, with in- 
clination ranging from 90 to 15 degrees, the resonant island ex- 
tends to large values of the rms of about 7 m/s, which reinforces 
the confidence that the present system is in a 3:1 resonant state. 

The present dynamical analysis also indicates that this 3:1 
mean motion resonance is stable over Gyrs timescale. We have 
as well tested directly this stability by performing a numerical 
integration of the systems Si and S2 (Tab. [B over 5 Gy r using 
the symplectic integrator SABAC4 of [Lask^& Robutell (1200 ID 
with a step size of 0.02 years. Both solutions <Si, <S2 remain sta- 



4 



J. Laskar and A.C.M. Correia: HD 60532, a planetary system in a 3:1 mean motion resonance 




Table 4. Change of the mass factor (1/ sin /) and libration period 
Pq with inclination (/). 



Fig. 3. Possible location of additional planet in the HD 60532 
system. For 90 (top) and 20 degrees (bottom) inclination, the 
stability of a small mass particle in the HD 60532 system is 
analyzed, for various semi-major axis and eccentricity, and for 
K = lO-^m/s. The stable zones where additional planets could 
be found are the dark blue regions. 



ble over 5 Gyr, although their orbital elements present large vari- 
ations because of the strong gravitational interactions between 
the two planets. For the nominal solution <Si, the eccentricity 
of the inner planet oscillates within 0.115 < et, < 0.301, while 
for the outer planet, we have 0.0165 < ec < 0.143. The secu- 
lar variations in the orbital parameters are mostly driven by the 
rapid secular frequency gi, of period 27r/gi ^ 754 yr (Table[3]). 
These secular variations in the orbital elements occur much more 
rapidly than in our Solar System, which may enable them to be 
detected directly from observations. 

5. Additional constraints 

The sta bility analysis summ arized in Figured] shows, as in HD 
45364 ( Correia et al.ll2009h . a very good agreement between the 
3:1 resonant islands and the rms contour curves. We can thus 
assume that the dynamics of the two known planets is not much 
disturbed by the presence of an additional large planet close by. 

We have tested the possibility of an additional planet in the 
systems Si and S2 by including a third planet d, with varying 
semi major axis and eccentricity over a large range, and per- 
forming stability analysis (Fig. [3]). The test was done for a fixed 
K value ; first with an asteroid size object (K = lO-^m/s) (Fig©, 
and then for an Earth-size object (K = O.lm/s) without notice- 
able change in the dynamics. 

From this analysis, one can see (Fig. [3]), that for Si (i = 90°), 
the zone of instability extends from about 0.2 AU to about 4 AU, 
while for S2 (i = 20°), no planet can be stable from 0.1 AU 
to about 6 AU. The eventual discovery of an additional planet 
in this range would thus constrain very much the inclination of 
the system with respect to the plane of view, and thus the plan- 
etary masses. At / = 90°, the libration period of the 3:1 reso- 
nant argument is 39.93 yr, but this period depends strongly on 
the planetary masses, and thus on the inclination / of the sys- 
tem (Tab.©. We have computed the libration period Pq for the 
best fit obtained at various inclinations, / = 90, 30, 20, 15, 10 de- 
grees. At / = 20°, we have Pq = 23.17 yr. As the amplitude of 
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Fig. 4. Radial velocity signature of the inclination (/). The diff'er- 
ence with respect to the nominal solution (at / = 90°), are plotted 
overtime, for / = 30° (red), 20° (green), 15° (blue), and 10° (vi- 
olet). The available HARPS residual data are also plotted with 
their error bars. 



the libration is large, we can expect that over the next decade, 
this resonant libration period will be constraint, allowing a good 
determination of the inclination and planetary masses. Indeed, if 
we extrapolate the diff'erences of radial velocity, for various in- 
clination hypothesis (/ = 30, 20, 15, 10 degrees), compared with 
the nominal solution Si, it appears clearly that the inclination 
of the system will be determined within a few years with the 
HARPS data (FigH. 



6. Discussion and conclusion 

We have re-analyzed the dynamics of the two planets system 
HD 60532, first reported in DES08. Contrary to the conclusions 
of this previous work, we believe that the global dynamical study 
presented here (Fig|4]) allows to conclude unambiguously that 
the system is in a 3:1 mean motion resonance. 

Moreover, due to the strong dynamical interaction present in 
this system, we could also fit the inclination which led to a most 
probable value of 20 degrees, while our dynamical analysis set 
up a stability constraints to 15° < / < 90°. The confidence in 
the presence of the 3: 1 resonance is enhanced by the fact that the 
system remains in resonance, for all the fits that we performed, 
from 90 to 15 degrees inclination. In the nominal solution {Si, 
Tabd]), the resonant angle 6 = At - 3 Ac -\- 2ojh is in libration 
around 180°, with a libration period of 39.93 years with an asso- 
ciated amplitude of 38? 14, although the main oscillation of this 
argument is driven by the secular term gi-g2 with an associated 
amplitude of 63? 65, the total amplitude of the libration being 
123? 16, due to additional harmonics. 

The planet-planet interactions in this system are large, due 
to the presence of this 3:1 resonance. The dependence of the 
libration period with the mass factor (Tab.© should should allow 
to determine the inclination of the orbital planes and thus the 
masses of both planets within the next decade (Fig|4l). 
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